Black Hole Production by Cosmic Rays 
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Ultra-high energy cosmic rays create black holes in scenarios with extra dimensions and TeV- 
scale gravity. In particular, cosmic neutrinos will produce black holes deep in the atmosphere, 
initiating quasi-horizontal showers far above the standard model rate. At the Auger Observatory, 
hundreds of black hole events may be observed, providing evidence for extra dimensions and the 
first opportunity for experimental study of microscopic black holes. If no black holes are found, the 
fundamental Planck scale must be above 2 TeV for any number of extra dimensions. 
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Black holes are among the most captivating and inac- 
cessible phenomena in physics. In principle, tiny black 
holes can be produced in particle collisions with center- 
of-mass energies above the Planck scale M*, where they 
should be well described semiclassically and thermody- 
namically jl[]. However, in conventional four-dimensional 
theories, M* ~ 10 19 GeV. Given currently accessible 
energies < 1 TeV, the study of such black holes is far 
beyond the realm of experimental particle physics. 

In models with extra dimensions, however, the funda- 
mental Planck scale may be much lower. If this is the 
case, black hole production and evaporation might be 
observed in particle collisions H, || . Beginning later this 
decade, for example, the Large Hadron Collider (LHC) 
at CERN will begin operation with parton center-of-mass 
energies of several TeV. Assuming M* of order 1 TeV, the 
authors of Refs. ^, [j| have noted the distinctive charac- 
teristics of black hole production and find event rates as 
large as 10 8 per year. 

As high as it is, the energy range probed by the LHC 
is modest compared to that of ultra-high energy cosmic 
rays, which have been observed to interact in the Earth's 
atmosphere with center-of-mass energies in excess of 100 
TeV. As we will see, cosmic neutrinos with energies above 
10 6 GeV are particularly effective sources of black holes, 
with production cross sections as much as two or more 
orders of magnitude above standard model (SM) predic- 
tions. These black holes decay rapidly, initiating spec- 
tacular quasi-horizontal showers deep in the atmosphere. 

Observation of such showers at the rates we predict 
would be a strong indication of new TeV-scale physics. 
In the SM, while Earth-skimming ultra-high energy neu- 
trinos may be observed with reasonable rates by fluores- 
cence detectors |6| and ground arrays , those that pass 
only through the atmosphere are extremely difficult to 
detect. Even at the large Pierre Auger Observatory, SM 
interactions are expected to produce only a fraction of an 
event per year |8[ ||, |l(|. In contrast, with conservative 
neutrino flux estimates, we find that Auger could detect 
hundreds of black holes before the LHC begins opera- 
tion, providing evidence for TeV-scale gravity and extra 



dimensions and making possible the experimental study 
of black holes in the late stages of Hawking evaporation. 

Low-scale gravity may be realized if the conventional 
four spacetime dimensions are supplemented by n addi- 
tional spatial dimensions. SM matter and gauge fields 
are typically assumed to be confined to the four dimen- 
sions of our world. However, gravity propagates in the 
full (4+n)-dimensional space with Einstein action 
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where M* is the fundamental Planck scale. If the 
?i-dimensional space is flat and compact with volume 
V n (TTJ , the observed gravitational strength is reproduced 
provided Ml +n V n « (1.2 x 10 19 GeV) 2 . For large V n , 
M* near the TeV scale is possible: current bounds are 
M* > 40 TeV for n = 2, M* > 400 GeV for n = 4, and 
M* > 300 GeV for n = 6 Jl|. For n = 1, M* - TeV is 
excluded in this context. However, in alternative scenar- 
ios with warped metrics fill , gravity may become strong 
below 1 TeV, even for n = 1 p4|] . Our discussion will not 
depend on the details of these scenarios, as long as the 
black holes produced are smaller than the compactifica- 
tion radii (curvature scales) in flat (warped) scenarios, 
and so well approximated by (4+n)-dimensional asymp- 
totically Minkowskian solutions. 

In theories with TeV-scale gravity, many effects may 
alter cosmic ray interactions with observable conse- 
quences ]i~5[ ]. However, in contrast to all other processes, 
our understanding of black hole properties is expected to 
be qualitatively sound for center-of-mass energies above 
M* and increasingly valid the farther above M* one 
probes. Black holes, then, may provide a uniquely reli- 
able probe of extra-dimensional effects above M* , where 
the full range of cosmic ray energies may be exploited. 

The Schwarzschild radius for a (4+n)-dimensional, 
neutral, non-spinning black hole with mass Mbh is JTH] 
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Black hole formation is expected when partons i and j 
with center-of-mass energy pass within a distance 
r s (s), suggesting a geometrical cross section of order 
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We will take this as an adequate approximation and as- 
sume that a black hole of mass -Mbh = V$ is formed. 
(Numerical analysis of classical head-on collisions in four 
dimensions finds Mbh ~ 0.8 VI Q.) The suppression 
factor of Ref. 24 has been disputed ^5|; we have not 
included it here. The neutrino-nucleon scattering cross 
section is then 



„(WV^BH) = W' 



dx&i(xs) fi(x,Q) , (4) 
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where s = 2ttinE u , the sum is over all partons in the nu- 
cleon, the fi are parton distribution functions (pdfs) , and 
M-gg 11 is the minimal black hole mass for which Eq. (|^) 
is expected to be valid. We set momentum transfer 
Q = min{A/BH, 10 TeVj, where the upper limit is from 
the CTEQ5M1 pdfs Jl§; <j{vN -> BH) is insensitive to 
the details of this choice. For the conservative fluxes 
considered below, our results arc also rather insensitive 
to x < 10~ 5 . For concreteness, however, we extrapolate 



to x < 10 assuming fi(x, Q) 



-[i+Ai(Q)] 



Finally, we 



choose -Mg™ = M*. The relatively mild dependence on 
Mg™ is discussed below. 

Cross sections for black hole production by cosmic neu- 
trinos are given in Fig. |l|. The SM cross section for 
vN — > IX is included for comparison. In contrast to 
the SM process, black hole production is not suppressed 
by perturbative couplings and is enhanced by the sum 
over all partons, particularly the gluon. In addition, 
while the SM cross section grows rapidly with E v , as 
is well known, the black hole cross section grows even 
more rapidly: for large n, it has the asymptotic behavior 
a cx E^ {w TeV) « E° A5 . As a result of these effects, 
black hole production may exceed deep inelastic scatter- 
ing rates by two or more orders of magnitude. 

Although greatly reduced by black hole production, 
neutrino interaction lengths L = 1.7 x 10 7 kmwe (pb/cr) 
are still far larger than the Earth's atmospheric depth, 
which is only 0.36 kmwe even when traversed horizon- 
tally. Neutrinos therefore produce black holes uniformly 
at all atmospheric depths. As a result, the most promis- 
ing signal of black hole creation by cosmic rays is quasi- 
horizontal showers initiated by neutrinos deep in the at- 
mosphere. At these angles, the likelihood of interaction 
is maximized and the background from hadronic cosmic 
rays is eliminated, since these shower high in the atmo- 
sphere. The number of black holes detected is, then, 



N = 



dE v N A —o{E v )A{E v )T, (5) 
ahj,. 



io 6 



io 4 



10 £ 




10 ! 



10° 10 1 

E„ (GeV) 



io 1 



FIG. 1: Cross sections a(yN -> BH) for M» = M^ 1 = 
1 TeV and n — 1, ... ,7 from above. (The last four curves are 
virtually indistinguishable.) The dotted curve is for the SM 
process vN — > IX. 



where A(E V ) is a given observatory's acceptance for 
quasi-horizontal showers in cm 3 water equivalent steradi- 
ans (cm 3 we sr), Na = 6.02 2 x 10 23 is Avogadro's number, 
d&/dE u is the source flux of neutrinos, and T is the run- 
ning time of the detector. 

There are many possible sources of ultra-high energy 
neutrinos. Here we conservatively consider only the 
'guaranteed' flux of Greisen neutrinos produced by in- 
teractions of the observed ultra-high energy cosmic rays 
with the cosmic microwave background |19(| . This flux is 
subject to uncertainties; we adopt the results of Ref. [po| , 
shown in Fig. 0. The flux estimates of Refs. |y1 produce 
similar event rates, while the strong source evolution case 
of Ref. [p2j enhances the results below by over an order 
of magnitude. New physics might also increase the neu- 
trino flux. In particular, many proposed explanations 
of cosmic rays with energies above the Greisen-Zatsepin- 
Kuz'min cutoff jl9|, would boost these event rates by 
several orders of magnitude. 

Quasi-horizontal showers may be observed by air 
shower ground arrays or air fluorescence detectors. The 
largest near-future cosmic ray experiment is the Auger 
Observatory, a hybrid detector consisting of two sites, 
each with surface area 3000 km 2 . Construction of the 
southern site is in progress, with a counterpart planned 
in the northern hemisphere. Auger acceptances for 
deeply penetrating air showers have been studied in 
Refs. 1 
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Black holes decay thermally, ac- 
cording to the number of degrees of freedom available, 
and so their decays are mainly hadronic |j| . We there- 
fore consider the hadronic shower acceptance for ground 
arrays, including 'partially contained' showers ||. For 
fluorescence, we use the results of Ref. Q for showers 
with zenith angles above 60° initiated at depths greater 
than 1250 cmwe. These acceptances are given in Fig. |^. 
A duty cycle of 10% has been included for fluorescence, 



E„ (GeV) 

FIG. 2: Neutrino flux from Greisen photoproduction 
(solid) H], and ground array (dashed) II and fluores- 
cence (dotted) |^] acceptances of one Auger site for quasi- 
horizontal hadronic showers. For fluorescence detection, a 
duty cycle of 10% has been included. 



where observations are limited to cloudless, moonless 
nights. At E u ~ 10 10 GeV where the Greisen flux peaks, 
the ground array is more sensitive, and so we focus on 
ground array rates below. Note, however, that future 
detectors, such as Telescope Array and the space-based 
OWL and EUSO, will improve this fluorescence accep- 
tance by one to three orders of magnitude [fj6f . 

Given the cross sections of Fig. [j] and the flux and 
acceptances of Fig. p| the number of events is determined 
by Eq. ([5]). The results for ground arrays are given in 
Fig. Tens to hundreds of events are possible for M* rs 

1 TeV. Tens of black holes may also be detected by 
fluorescence. For larger M*, — > BH) falls rapidly 
as A/* ( 4 + 2n )/( 1 +™)^ Nevertheless, requiring 3 events for 
discovery, black hole production probes Planck scales as 
high as 3 TeV for n = 1, and 2 TeV for all n. If no events 
are seen, barring a neutrino flux significantly below our 
conservative estimate, a stringent lower bound of M* > 

2 TeV may be set for all n. 

The results of Fig. | are for M™ n = M*. While the 
semiclassical approximation is invalid for Mbh ~ 
this is a calculational, not physical, limitation and does 
not imply that production of black holes or similar states 
in this mass range is suppressed p5| ; in fact, it may just 
as well be enhanced. Nevertheless, it is comforting to 
know that our results are not strongly sensitive to this 
assumption. In Fig. ||, the dependence on Mg 1 ^ 11 is shown. 
For Mg™ = 5M«, event rates are reduced by factors of 
2 for n = 1 and 4 for large n. While these reductions are 
substantial, they are extremely mild relative to the case 
at colliders. At the LHC, the requirement A/bh > 5M* 
suppresses event rates by factors of a hundred or more |j) . 
For cosmic rays, while the black hole mass distribution 
is still peaked at low masses as a result of enhancements 
from pdfs at low x, the reduction is far more modest. 
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FIG. 3: The number of black holes detected by the ground 
array in 5 Auger site-years as a function of M* = Mbh" and 
the number of extra dimensions n. 
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FIG. 4: The number of black holes detected by the ground 
array in 5 Auger site-years as a function of M^ 1 for M* 
1 TeV and n = 1, . . . , 7 from above. 

If an anomalously large quasi-horizontal shower rate 
is found, it may be identified as due to black hole pro- 
duction in several ways. First, although a large rate 
may be attributed to either an enhanced flux or an en- 
hanced black hole cross section, these possibilities may 
be distinguished by searches for Earth-skimming neutri- 
nos H 0. While an enhanced flux increases these rates, 
a large black hole cross section will suppress them, since 
the hadronic decay products of black hole evaporation 
will not escape the Earth's crust. 

Second, showers from black hole production have dis- 
tinctive characteristics. In the SM, typical hadronic 
showers, as initiated by nucleons or nuclei, occur high 
in the atmosphere. Deep atmospheric showers arise only 
from vN — > £X, resulting in a hadronic shower initi- 
ated by the struck quark, possibly accompanied by an 
electromagnetic shower carrying most of the incident 
energy, depending on the neutrino flavor. Black hole 
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events are markedly different. The black hole rest life- 
time is r ~ (l/M*)(M B n/M*)( 3+n V ( - 1+n \ Since M" 1 ~ 
TeV -1 ~ 10~ 27 s, even the largest black holes produced 
evaporate effectively instantaneously. In contrast to SM 
showers, however, black hole showers have small elec- 
tromagnetic components, and the average multiplicity in 
black hole decays is ||) 
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where Tjj is the Hawking temperature. Large mass black 
holes therefore decay to large numbers of quarks and glu- 
ons, and black hole showers will appear more nucleus-like 
than SM events, with the discrepancy growing with black 
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in s everal wayo pi f r — X max , the atmoophoric depth at 

which the number of particles in a shower reaches its 
maximum, is significantly lower for nuclei, and shower- 
to-shower fluctuations in A max and the number of elec- 
trons are also smaller. Black holes and SM events may 
therefore be distinguished based on shower characteris- 
tics, at least on a statistical basis. Note from Fig. [| that 
a fairly smooth distribution of black hole masses is ex- 
pected. If large numbers of black holes are found, the 



correlations of shower energy with AIbh and A max with 



(N) will also allow tests of Hawking evaporation and pos- 
sibly even measurements of n and M* through Eq. (0). 

Before closing, we comment on the possible relevance 
of black hole production to the GZK paradox. As noted 
above, the cross sections of Fig. [l] may be enhanced, es- 
pecially for Mbh ~ M*, where the behavior of black 
holes and related objects is very poorly understood. In 
addition, if effectively four-dimensional black holes are 
produced, as may be possible in warped scenarios with 
small curvature scales, we find cross sections of 10 mb for 
E v ~ 10 12 GeV. If these or other enhancements are large 
enough to bring the cross sections to the 100 mb level, 
cosmic neutrinos will be primaries immune to GZK-type 
cutoffs that produce hadronic showers high in the atmo- 
sphere, providing a viable resolution to the GZK puzzle. 
While the required enhancement is large and speculative, 
the qualitative merits of black hole production as a so- 
lution to the GZK paradox are suggestive and deserve 
further study. 

To summarize, in TeV-scale gravity scenarios, ultra- 
high energy cosmic neutrinos will produce black holes 
in the Earth's atmosphere, leading to anomalously large 
rates for quasi-horizontal hadronic showers. If the LHC 
is to be a black hole factory, at least tens to hundreds of 
black holes will be detected at the Auger Observatory be- 
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